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When two sheets of graphene stack in a twisted bilayer graphene (tBLG) configuration, the result¬ 
ing constrained overlap between interplanar 2 p orbitals produce angle-tunable electronic absorption 
resonances. Using a novel combination of multiphoton transient absorption (TA) microscopy and 
TEM, we resolve the resonant electronic structure, and ensuing electronic relaxation inside single 
tBLG domains. Strikingly, we find that the transient electronic population in resonantly excited 
tBLG domains is enhanced many fold, forming a major electronic relaxation bottleneck. 2-photon 
TA microscopy shows this bottleneck effect originates from a strongly bound, dark exciton state 
lying ~0.37 eV below the 1-photon absorption resonance. This stable coexistence of strongly bound 
excitons alongside free-electron continuum states has not been previously observed in a metallic, 2D 
material. 


Photoexcited electrons in graphene relax energeti¬ 
cally far faster than the e-h separation timescale, mak¬ 
ing many electronic and optoelectronic applications 
prohibitive. m While similar fast, picosecond relax¬ 
ation timescales are also observed in Bernal stacked bi¬ 
layer graphene (FBLG). slower relaxation might be pos¬ 
sible in twisted bilayer graphene (tBLG). In tBLG, an 
off-axis interlayer twist angle ( 0 ) gives rise to band anti¬ 
crossings and van Hove singularities (vHs, Fig lb).j6HS] 
Near such r>Hs, previous studies show that optical ab¬ 
sorption increases by ^20% and is peaked at an energy, 
EY [61171I9HT6] This absorption resonance peak increases 
monotonically with 6 (see supplementary video). [IT] To 
date however, the properties of photoexcited electrons 
in vHs remain unexplored beyond the Raman and linear 
absorption characterization. Here, we apply space, time, 
and energy-resolved 1-photon (1-ph) and 2-photon (2-ph) 
transient absorption (TA) microscopy to both spectrally 
map the excited state electronic-structure of tBLG and 
image the ensuing electronic dynamics. 



FIG. 1. (a) tBLG free-electron interlayer band-structure, 

(b) Cross-sectional view shows interlayer rHs resonances, X 13 
and X 24 , between band anti-crossing regions, (c) Alterna¬ 
tively, the degenerate X± 3 and X 2 4 states may rehybridize, 
giving a 1-ph state above, and a 2-ph allowed exciton state 
below. TA (arrows) interrogates the electronic population 
(circles). 


The single-particle band structure for tBLG can be un¬ 
derstood by superimposing two graphene Brillouin zones, 
rotated by a twist angle 0 , as shown in Fig. la. mum 
The vertical line cutting through the two Dirac points of 
the graphene layers (Fig. lb) shows the band anticross¬ 
ing near the degeneracy with an energy splitting (A), 
and four possible vHs transitions between the graphene 
sub-bands labeled 1 through 4 . These optical transitions 
experience a large joint density of states between the 
valance bands (1 & 2 ) and the conduction bands (3 & 
4 ), but only 1 —>>3 (denoted X13) and 2— >4 (X 24 ) transi¬ 
tions are allowed due to selection rules. [ 111120 ] 

Outwardly, the X13 and X24 transitions shown in Fig. 
lb are degenerate vHs similar to graphene’s M -point 
saddle-point exciton. nami In this case, Coulombic at¬ 
traction between e-h pairs would augment both X13 and 
X24 transition energies and produce an asymmetric, Fano 
optical lineshape.[22, 23 ] Since such unbound Fano ex¬ 
citons couple to continuum states of graphene. 22 this 
model predicts the 1-ph TA response from tBLG will de¬ 
cay quickly, with an amplitude and rate similar to single¬ 
layer graphene (dotted arrow in Fig. lb). 

Alternatively, previous studies suggest that inclusion 
of bound-exciton effects are necessary to simulate the 
nearly gaussian tBLG optical absorption lineshape.[ 17 , 
[ 24 ] While one can consider unbound excitonic states for 
X13 and X24 independently, such a picture is incomplete 
because the two states occur at the same energy and 
momentum, a direct result of the electron-hole symme¬ 
try at the r>Hs in tBLG. A more complete description 
was given in recent work reported by Liang et al. that 
predicts formation of stable, strongly bound ( Eb ^ 0.5 
eV) excitons. [ 24 ] These first-principles calculations sug¬ 
gest that interlayer excited states in tBLG are better de¬ 
scribed by renormalized symmetric (X$=X 13+Y24) and 
anti-symmetric (Xa=X 13-Y24) excitonic states. [22 In 
this model, illustrated in Fig. lc, X$ corresponds to 
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the optical tBLG resonance at Eq , and is an unstable 
exciton. min] Conversely, Xa is only 2-ph accessible and 
is calculated to be a strongly-bound, localized excitonic 
state. [24 The remarkable stability predicted for the Xa 
state results from the deconstructive coherence between 
the two degenerate Fano resonances rigorously cancel¬ 
ing coupling with graphene continuum states. Such a 
state is termed a ’ghost Fano’ resonance. While similar 
phenomena have been weakly observed in quantum dot 
and carbon nanotube systems, such strongly bound ex¬ 
citon states have never been observed in a 2D metallic 
system. [Ml [26 If such ghost Fano excitons are present in 
tBLG, weak exciton-continuum coupling is expected to 
enhance the local electronic population, giving longer re¬ 
laxation dynamics for the 1-ph (X$) and 2-ph (Xa) TA 
response. 

In this work, we obtain the TA spectra and dynamics 
of single tBLG domains, and map out the different 1-ph 
and 2-ph electronic transitions predicted by the contrast¬ 
ing vHs and strongly-bound exciton models in Fig. 1. We 
further correlate ultrafast TA microscopy with the precise 
local atomic stacking and grain boundaries, by employ¬ 
ing darkfield TEM to definitively assign a twist angle to 
the absorption resonance, Ea.\ 27] Our experimental TA 
microscopy approach is outlined in Fig. 2a. The 1-ph TA 
map in Fig. 2a shows a prominent patch of 6.8° oriented 
tBLG that is surrounded by non-twisted CVD graphene 
on a silicon nitride membrane substrate. This map was 
obtained by raster scanning a diffraction-limited pump 
and probe pulse pair over the graphene. We tuned our 
140 fs pump pulse to be resonant with the 6.8° domain 
at E purnp ~Eo~ 1.3 eV. After a delay time t, we detect 
the differential TA (A R(t) oc A cr(t)) of a collinear probe 
pulse and construct time-dependent TA maps point-wise. 
Using probe energies (E pr = 0.8 eV) well below the reso¬ 
nance Eq, in Fig. 2a graphene gave an interband de¬ 
creased absorption response (i.e. Pauli blocking of probe 
beam) everywhere at all time delays. [2] 

Our TA maps can be interpreted as ’movie frames’ that 
closely approximate the relative photoexcited electronic 
population at a particular probe energy and time-delay 
(see supplementary video 2). The 6.8° tBLG region la¬ 
beled in Fig. 2a has a ^two-fold stronger TA Pauli block¬ 
ing response than the adjacent 0° stacked regions. How¬ 
ever, the corresponding linear absorption map in Fig. 2c 
only shows a ^20% resonant enhancement. To account 
for this discrepancy, electrons in interlayer tBLG avoided 
crossing regions must relax much slower than the sur¬ 
rounding non-twisted graphene bilayers, suggesting an 
intrinsic electronic relaxation bottleneck. 

In Figure 2b, we repeat the 1-ph measurement but 
instead resonantly probe the electronic population at 
Eq (E pump = 1.33 eV, E pr = 1.26 eV probe). Compared 
against the corresponding linear absorption map in Fig. 
2c (inset), the TA maps differ in both sign and absolute 
amplitude. Strikingly, only the 6.8° tBLG domain gives 
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FIG. 2. tBLG electronic relaxation bottleneck (a) Ul- 

trafast scanning TA microscopy map of multilayer graphene. 
The TA Pauli blocking response is >2x enhanced for the 
6.8° tBLG domain, (b) On-resonance TA maps at E purnp ~ 
E pr = Eq show a TA response localized to the 6.8° domain 
is present for t>100 ps . The surrounding graphene gives 
only a weak transient intraband response (opposite sign), (c) 
TA relaxation kinetics of the 6BLG and tBLG regions labeled 
(vs. SWNT Eh state, gray). Corresponding linear absorp¬ 
tion map at 1.3 eV shows tBLG (red) is only ~20% stronger 
than in FBLG (inset). 


a strong TA Pauli blocking response. Meanwhile, the 
surrounding graphene in Fig. 2b gives a weak, short¬ 
lived graphene intraband TA response signified by its 
opposite sign. This suggests that interlayer tBLG elec¬ 
trons are decoupled from the intraband transient re¬ 
sponse that dominates the TA map everywhere else in 
Fig. 2b. Surprisingly, the subsequent TA movie frames 
show excited carriers are present even >100 ps after ini¬ 
tial excitation. Both observations definitively show that 
interlayer electrons excited at Eq experience a major bot¬ 
tleneck restricting electronic relaxation. Such a strong 
and long-lived electronic signal in tBLG disagrees with 
the continuum Fano resonance model (Fig. lb), but 
can be explained by an excitonic model (Fig. lc) where 
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FIG. 3. 1 photon vs. 2 photon absorption of ~6.5° and ~8°(dotted outlines) tBLG domains, (a) Competing models; 
(z.)bound exciton model, (ii .)continuum model (b) TA spectrum (red) vs. 2-ph transition pump energy. 1-ph interlayer 
absorption spectrum (crtBLG-crbBLGi green) for 8° tBLG. The lowest 2-ph peak fits best to a gaussian lineshape centered at 
5 =0.37 eV below Eq. There is also a Fano resonance lineshape feature at A = 0.33 eV above. ( inset , map of tBLG absorption 
resonance vs. twist angle) (c) 1-photon TA map, at E purnp = 1.3 eV shows a strong electronic bleach only from the resonantly 
excited 6.5° domains, (d) Conversely, a 2-ph TA map at E vurnp — 0.6 eV shows a ground-state bleach only from the 8° domains. 
Combined, these maps demonstrate that the (Xs) state is 1-ph allowed and (Xa) state is only 2-ph allowed. 


weak exciton-continuum coupling allows for stable exci¬ 
ton formation. m 

To isolate the relaxation rates intrinsic to interlayer 
tBLG electrons excited at Eq , we plot tBLG-5BLG (red) 
in Fig. 2c by subtracting the much weaker (and op¬ 
posite signed) intralayer electronic TA response (blue). 
A similar approach has been previously used to decou¬ 
ple linear absorption spectra, as cr tB LG~^bBLG • m ng 
A least-squares deconvolution exponential fit of the ki¬ 
netic decay requires only a biexponential fit that decays 
with lifetimes of 1.4T0.1 ps and 66±4 ps. These lifetime 
components are remarkably long for any electronic state 
within a metallic system. By repeating our linear and TA 
measurements at low temperatures, we further found the 
interlayer electronic response appears largely invariant 
to both lattice temperature (5-295 K) and the substrate 
used (see supplemental materials). Furthermore, the TA 
response did not shift sign as the probe wavelength was 
scanned through the resonance, Eq. Combined, these ob¬ 
servations suggest that laser-induced heating effects do 
not contribute appreciably to the overall large TA signal 
response, and so the TA signal is predominately elec¬ 
tronic in origin. 


If carrier in tBLG are unbound excitons, the TA (red) 
in Fig. 2c must relax at a rate similar to 6BLG ( navy ), 
providing that phonons with E> A (dotted arrow in Fig. 
lb) are available to scatter carriers through the anti¬ 
crossing gap (A) illustrated in Fig. lb. [2] Comparison 
of the short-time tBLG kinetics against 6BLG in Fig. 
2c (red) reveals the absence of the dominant fast sub-ps 
electron relaxation components associated with graphene 
electron thermalization and optic phonon emission. m 
Remarkably, the shortest interlayer tBLG lifetime is 1.4 
ps, which is similar to graphene’s rate-limiting relax¬ 
ation rate that is often associated with disorder-assisted 
or supercollision relaxation. [2j [28] The absence of the 
sub-ps relaxation processes, and the emergence of this 
long ^66 ps decay in tBLG relaxation kinetics, suggests 
that some electrons are decoupled from graphene’s con¬ 
tinuum states, as predicted by the strongly bound exciton 
model. [24] 

The unexpected TA bottleneck we observe in tBLG 
may be compared against semiconducting single-walled 
carbon nanotubes (SWCNTs), a similar carbon system 
with constrained 2 p orbital interactions. It is estab¬ 
lished that SWCNTs have 1-ph and 2-ph excitonic states 
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resulting from (chiral) angle dependent overlapping 2 p 
orbitals. |29h3T] Fig. 2c {gray) directly compares the 
Ei i exciton relaxation rate of (6,5) chirality SWCNTs 
against £BLG {red). While the short time behavior dif¬ 
fers greatly, Fig. 2c shows the longer components of both 
traces decay at a similar rate, suggesting that the dy¬ 
namic phonon environment causing En exciton relax¬ 
ation in SWCNTs might be of a similar nature to the 
interlayer exciton-phonon interactions causing exciton re¬ 
laxation in £BLG. 

We can better distinguish between competing uHs and 
bound exciton models outlined in Fig. 3a, by exploit¬ 
ing the 2-ph selection rules required for the predicted 
dark £BLG exciton state, Xa- pnnn] To search for pos¬ 
sible dark state transitions, we used a different sample 
of CVD bilayer graphene. The linear absorption map 
shown in Fig. 3b {inset) reveals a series of £BLG domains 
with twist angles of either ^6.5° {yellow) or ^8° (red), 
corresponding to Eq ^1.25 eV and 1.55 eV respectively. 
Figure 3b plots the 2-ph TA spectrum obtained using 
IR pump energies ranging from E purnp =0.49 to 1.15 eV, 
and a 8° resonant probe at Eq~E pv = 0.56 eV. We observe 
two clear TA peaks centered at 1.18 eV and 1.82 eV that 
originate from a resonantly enhanced 2-ph absorption. 
Specifically, as illustrated in Fig. 3ai, we observe these 
dark states through resonant 2-ph enhanced Pauli block¬ 
ing of the depleted ground state. Moreover, our ability 
to probe electronic population of optically dark state re¬ 
quires that Xs and Xa states share a common ground 
state; an inherent feature of a bound exciton model. m 

The lowest peak in Fig. 3b indicates that enhanced 
2-ph absorption took place via a discrete, low-lying tran¬ 
sition centered at 1.18 eV. Comparing the 2-ph peak 
against the 1-ph absorption resonance at Eq= 1.55 eV 
(green, Fig. 3b), we readily obtain the energy-state split¬ 
ting parameters of S = 0.37 eV and A = 0.33 eV. This 
0.37 eV energy splitting closely matches the theoretically 
predicted S =0.4-0.5 eV, state splitting calculated for 21° 
tBLG.[24] Such a large bright-dark state energy split¬ 
ting is much greater than the analogous state splitting 
in SWCNTs, and explains why photoluminescence has 
not yet been observed from resonantly excited £BLG do¬ 
mains. 

To completely map the selection rules associated with 
£BLG electronic transitions, we compare the 1- and 2-ph 
TA microscopy response of 6.5° and 8° oriented domains 
in Fig. 3c-d. Fig. 3c maps-out the 1-ph TA response for 
Epump— 1-3 eV, E pr = 1.2 eV at £=0.3 ps. Despite the Xa 
state being 1-ph resonant with the ^8° £BLG (dotted red 
outlines), we observed only a weak intraband response 
as was seen for 6BLG regions previously (Fig. 2b). This 
confirms that the Xa transition is not 1-ph accessible. 
In contrast, the 6.5° £BLG domains give a strong Pauli 
blocking response because E pump is resonant with Xs- 

2-ph resonant transitions of single £BLG domains are 
imaged in Fig. 3c, by tuning our pump pulse energy to 


roughly half the predicted 8° Xa state energy(see Fig. 
3a), or E pump = 0.6 eV. Comparison of the TA maps in 
Fig. 3c against 3d show all of the 8° £BLG domain exci¬ 
tations that were forbidden under 1-ph excitation condi¬ 
tions are now allowed for a 2-ph excitation. Conversely, 
all the 6.5° £BLG domain excitations that were observed 
under 1-ph excitation conditions now appear dark (inac¬ 
cessible) under two-photon excitation. Using state parity, 
we assign the two-photon accessible dark states in Fig. 
3d to electronic carriers populating the Xa state of 8° 
£BLG. Together, Fig. 3c and 3d show that the Xs bright 
state is two-photon forbidden, and the dark Xa state 
is only two-photon allowed. These strongly enforced se¬ 
lection rules follow the parity expectations of a roughly 
hydrogenic-like, strongly-bound exciton model advocated 
by recent first-principle simulations. m 

The 2-ph spectral peak centered at 1.82 eV in Fig. 
3b has not been previously predicted or observed. This 
peak has a broader, asymmetric shape, that fits better 
to a Fano lineshape expected from the unbound exci¬ 
ton model. [as mm In contrast, the other 1-ph and 2-ph 
peaks in Fig. 3b fit best to a gaussian lineshape, a com¬ 
mon characteristic of bound excitonic transitions. Ac¬ 
cordingly, we infer that the two-photon absorption near 
1.82 eV is best assigned to an unbound state transition 
labeled X $+A in Fig. 3ai. Conversely, the two-photon 
absorption resonance lying S = 0.37 eV below is best 
characterized as the Xa bound exciton or ghost Fano 
resonance peak predicted by Liang et al. [24] as sup¬ 
ported by (i.) its asymmetric energy spacing (i.e. S vs. 
A), (ii.) gaussian lineshape, (iii.) long electronic lifetime 
and (iv.) parity enforced two-photon selection rules for 
the Xa and Xs transitions. 

Using both explicit calculations based on the Bethe- 
Salpeter equation and effective low-angle continuum 
model, Liang et al. predicted radically different electron¬ 
ics properties emerge for both the theorized Xs and Xa 
exciton states. [24] Specifically, the symmetric Xs state 
was found to have delocalized wavefunctions and a negli¬ 
gible binding energy. Conversely, the antisymmetric state 
Xa is predicted to be optically dark, insensitive to e — h 
charge screening, and strongly bound. [24] While certain 
phonons can scatter bound excitons into the lower-lying 
continuum states, the exciton-continuum coupling for the 
Xa is predicted to be vanishingly small, and roughly in¬ 
tensive to charge screening effects. m Accordingly, both 
theory and our TA microscopy now support that fast ex¬ 
citon dissociation becomes unfavorable in the Xa state 
of £BLG, enabling stable and meta-stable bound exciton 
states to form. 

Lastly, we compare the 1-ph Xs {pink) vs. 2-ph Xa 
{ black) electron relaxation kinetics measured for 8° £BLG 
in Fig. 4. We find that the normalized TA relaxation 
kinetics for one-photon and two-photon resonant excita¬ 
tions are nearly identical. These matching kinetics indi¬ 
cate that both signals originate from the same depleted 
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FIG. 4. Resonant vs. non-resonant electronic re¬ 
laxation for 8° tBLG. The vertical pump (solid arrows ) and 
probe (dashed arrows) combinations labeled on the absorp¬ 
tion spectrum (inset) of 8° tBLG , correspond by colour to 
the normalized 1- or 2-ph. relaxation kinetics plotted. 


common ground state, and that the electrons are im¬ 
pulsively transferred from the bright (Xs) to the dark 
(Xa) state as illustrated in Fig. 3ai. Such fast kinetic 
Xs — >Xa relaxation is consistent with theory showing 
that Xs is an unstable exciton state. [24] As a control, in 
Fig. 4 (gray) we show that the relaxation kinetics are 
impulsive when both E purnp and E pr > Eq , indicating 
that only free electron states are probed above resonance. 
Lastly we consider the case of E pr < Eq < E purnp , and 
find the long-decay components are consistently absent, 
suggesting resonant optical excitation may be required 
to form a long-lived stable exciton. Nonetheless, the off- 
resonance tBLG kinetic relaxation (green) is still signif¬ 
icantly enhanced in amplitude and lifetime compared to 
the bBLG TA response (Fig. 2c). This suggests that an 
electronic relaxation bottleneck effect is still present even 
when tBLG is optically excited above resonance. 

By definitively isolating the interlayer electronic dy¬ 
namics of 1- and 2-ph resonant optical transitions in 
tBLG, we have uncovered a fine-structure of bound (Xa) 
and unbound (Xs, X^-fA) exciton states that agrees well 
with recent simulations. 21j Specifically, we employed a 
novel form of diffraction-limited TA microscopy to ob¬ 
tain the intrinsic spectra and dynamics of single tBLG 
domains under a variety of resonant and non-resonant 
pump/probe combinations. In Fig. 2b, we show a TA- 
movie of electronic population that reveals the striking 
contrast between the bound-exciton carriers in the tBLG 
region and the free-electron population in the surround¬ 
ing graphene. These results suggest that the photoex- 


cited tBLG interlayer electrons are initially decoupled 
from scattering into graphene continuum sates, and ex¬ 
perience a significant electron relaxation bottleneck. In 
particular, resonantly excited carriers in tBLG give much 
stronger TA amplitude with longer relaxation kinetics for 
both the short (~2 ps) and long relaxation timescales 
(^70 ps). This bottleneck is best explained by the exis¬ 
tence of strongly-bound excitons in a ’ghost Fano’ state 
that we explicitly resolve using 2-ph TA microscopy. (24] 
Our results imply that tBLG may be a unique hybrid 
electronic material where free-electron metallic character 
can coexists alongside stable exciton states. The work 
further opens up possible new avenues for carrier extrac¬ 
tion that combine the high conductivity of metallic in- 
tralayer-electrons, with the enhanced electronic popula¬ 
tion that is now established for the interlayer electrons 
in tBLG. 

Methods : Multi-layer graphene was grown using low 
pressure CVD method on copper foil and transferred to 
silicon nitride grids (see supplemental materials).[ 34] Ar¬ 
eas containing low-angle tBLG were first identified using, 
a combination of hyperspectral absorption imaging tech¬ 
nique, and dark field TEM (DF-TEMb[TT] Final twist 
angle assignments of the bilayer patches were made by 
correlating the linear absorption and 1-ph TA peak spec¬ 
tral peaks energies. 

tBLG bright (Xs) and dark (Xa) states and their cor¬ 
responding electronic dynamics were measured using 1- 
and 2-ph confocal scanning TA microscopy. [35] Collinear 
pump-and probe pulses were obtained from two indepen¬ 
dently tunable outputs of an ultrafast system composed 
of Ti:Saph oscillator (Coherent Chameleon Ultra II, 80 
MHz, wavelength range 680-1080 nm) pumping an op¬ 
tical parametric oscillator (APE-Compact, wavelength 
range 1000-4000 nm). For one-photon TA measurements 
requiring pump and probe pulse doubly resonant with 
the bright (Xs) transition, a white-light supercontinuum 
probe was instead used. Cross-correlation of the pump 
and probe after the objective yielded a FWHM pulse du¬ 
ration of 142 fs. 

After a mechanical delay stage, both the pump, and 
the probe beams were aligned in a collinear geometry, 
raster-scanned by piezo-scanning mirror and coupled into 
a confocal scanning microscope via a 50X IR-region en¬ 
hanced, achromatic objective (NA= 0.65). One- and 
two photon transient absorption signals were detected 
by measuring the probe beam on with a TE cooled In- 
GaAs detector connected to a Zurich HF2LI lock-in am¬ 
plifier. The pump beam was modulated at either 0.25 
or 1 MHz using a AO-modulator (Gooch & Housego) 
to enable high-frequency lock-in detection of the differ¬ 
ential reflectivity. Appropriate optical filters were used 
in front of the detector to block the pump beam. The 
pump and probe spot sizes on the sample were deter¬ 
mined to ^1.5 /im, by fitting to a confocal scanning re¬ 
flection profile of deposited gold pads. The fluence of the 
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probe power was 5% of the pump fluence. Except where 
specified, all the measurements were done at 295 K. The 
probe power was fixed at 1 x 10 12 photons/cm 2 ) for 
the pump power dependence measurements. Microscope 
objective/transmission corrections curves were measured 
and rigorously taken into account for all the wavelengths, 
after each measurement. 
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